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We discuss the techniques to extract the electromagnetic A form factors in Lattice QCD. We 
evaluate these form factors using dynamical fermions with smallest pion mass of about 350 MeV. 
We pay particular attention to the extraction of the electric quadrupole form factor that signals 
a deformation of the A. The magnetic moment of the A is extrapolated using a chiral effective 
field theory. Using the form factors we evaluate the transverse density distributions in the infinite 
momentum frame showing deformation in the A. 
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1. Introduction 

Electromagnetic form factors probe the structure of hadrons, yielding information on their size, 
shape and magnetization. While the nucleon form factors and the N ^ ^ transition form factors 
have been studied quite thoroughly both experimentally and on the lattice much less has been 
done for the A form factors. Experiments are notoriously difficult due to the short mean life time of 
the A of only about 6 x lO^^'*^. Nevertheless the magnetic moments of the A+ and A++ [||, ^] 
have been measured. Lattice Quantum Chromodynamics (QCD) provides a well-defined frame- 
work to directly calculate the A form factors from the fundamental theory of strong interactions. 
A primary motivation for this work is to understand the role of deformation in hadron structure. 
For hadrons a classical non-relativistic description is not adequate and the definition of shape has 
to be refined [^. A suitable framework to discuss the shape of such systems is to define transverse 
charge distributions in the infinite momentum frame. Thus, a major achievement of this work is 
the development of lattice methods with sufficient precision to extract the electric quadrupole form 
factor and connect it in a rigorous manner to a well-defined charge density distribution so that the 
shape of the A can be discussed. In order to evaluate the A electromagnetic (e.m.) form factors to 
the required accuracy, we isolate the sub-dominant electric quadrupole form factor (FF) from the 
dominant electric and magnetic EEs. This is crucial since without the construction of an optimized 
source for the sequential propagator the electric quadrupole form factor can not be extracted to the 
desired precision. We note that this increases the computational cost since additional inversions 
are needed. These techniques were first tested in quenched QCD [Q]. Here we present results in 
the quenched approximation, using two degenerate flavors {Np = 2) of Wilson fermions and in a 
mixed action approach where domain wall valence quarks on a staggered A^^ = 2 + 1 sea are used. 



2. The transverse charge densities for a spin-3/2 particle 

We define a quark charge density for a spin-3/2 particle, such as the A(1232), in a state of 
definite light-cone helicity A, by the Eourier transform |0, ^], 

in the Breit frame, where An (G^) are helicity amplitudes and = —q^ = {Pf — PiY with pi and 
p f the initial and final momentum. The two independent quark charge densities for a spin-3/2 state 
of definite helicity are given by pf (^) and pf{b). Note that for a point-like particle, the 'natural' 

2 2 

values lead to As 3 (2^) = 1 [^], implying Pi (^) = 5^(b). The above charge densities provide us 
with two combinations of the four independent A EEs. To get information from the other EEs, we 
consider the charge densities in a spin-3/2 state with transverse spin. We denote this transverse 
polarization direction by S± = cos05ev + sin05ey, and the A spin projection along the direction of 
S± by s±. We can define the charge densities in a spin-3/2 state with transverse spin s± as : 

PT.,(b) = / ^e-'^^^''-^{P\^4,s^\J^0)\P\-^4,s^). (2.2) 



(271)2 2^+' > 2 ' ' 2 

By expressing the transverse spin basis in terms of the helicity basis for spin-3/2 and working 

3 
2 



out the Eourier transform in Eq. (2.2) for the two cases where sj_ = ^ and sj_ = ^ one obtains: 
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dQ 
2n 



Q 



■/o(e^)^ (a33 + 3Ai i) - sin(0fo-0s)7i(e^)^ , + 3A,_/ 



/o 1 

cos[2((/)fe - (^sViiQb) ^A3_ . +sin[3((/)^, - h)]h{Qb) -A,_^_ 



(2.3) 



and 



,A n.. l^°"dQ „r 1 



PriW = / ■/o(e^)7 3A3 3+A 
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sin(0/,-0s)/i(eZ^)^ f2^/3A3i -Ai_i 
4V 22 22 



V3 



COs[2((/)fo-(/)5)]/2(e^) ^A3_l 



-MXh-h)]h{Qb) -A3_3 



(2.4) 



where we defined the angle 0^ in the transverse plane as, b = ^(cos^fo^v + sin^fog^). One notices 



from Eqs. ( |2.3|j2.4t ) that the transverse charge densities display monopole, dipole, quadrupole, and 
octupole field patterns, which respectively are determined by the helicity form factors with zero, 
one, two, or three units of helicity flip between the initial and final A states. 

We can now evaluate the electric quadrupole moment corresponding to the transverse charge 
densities Pj^ . Choosing = Cx, the electric quadrupole moment can be defined as : 



(2.5) 



Qf^^e d'bibi-b'y)p^,^{b). 



From Eqs. (2.3 2.4) one obtains 



-ef = 7 {2 [Gmi (0) - 3eA] + [Ge2{0) + 3eA]} ( 
2 2 



(2.6) 



We may note that for a spin-3/2 particle without internal structure, for which the tree-level pre- 
dictions are Gmi(0) = 3eA and Ge2{^) = — 3eA, the quadrupole moment of the transverse charge 



densities vanishes. It is thus interesting to observe from Eq. ( ^^ ) that Qf^ is only sensitive to the 
anomalous parts of the spin-3/2 magnetic dipole and electric quadrupole moments, and vanishes 
for a particle without internal structure. The same observation was made for the case of a spin- 1 
particle in Ref. [^]. Furthermore, the factor 1/2 multiplying the curly brackets on the right hand 
side of Eq. ( |2.6| ) can be understood by relating the quadrupole moment of a 3-dimensional charge 
distribution (which we denote by Q^^) to the quadrupole moment of a 2-dimensional charge distri- 
bution (denoted by Q2d) both defined w.r.t. the spin axis. By taking the spin axis along the ;c-axis, 
the quadrupole moment for a 3-dimensional charge distribution pj^ is defined as : 

Q3d = J dxdydz {3x^ - r^) P3d{x,y,z),= J dxdydz [{x^ - /) + - z^)] P3d{x,y,z). (2.7) 

For a 3-dimensional charge distribution which is invariant under rotations around the axis of the 

^2 



spin, the two terms proportional to {x^ —3'^) and (x^ — z^) in Eq. (2^) give equal contributions 
yielding : ^3^/ = 2/ dxdydz {x^ —y'^)p^d{x,y,z)- Introducing the 2-dimensional charge density in 
the x3'-plane as : p2d{x,y) = J dzp3d{x,y,z), one immediately obtains the relation Q^d = '^Qid, 
with the quadrupole moment of the 2-dimensional charge density defined as : Q2d = J dxdy [x^ — 
y^) P2d{x,y)- Because Q^d is proportional to Ge2{0) in our case, we see that Q2d, which is half the 



value of Ge2{0), is consistent with Eq. (2.6) 
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We can also evaluate the electric octupole moment corresponding to the transverse charge 
densities Pj^^. Choosing = e^, the electric octupole moment can be defined as : 

0^^=ej cfbb^ sm{3^t)p^,,^ (b) = e J cfbby {3b^ - b]) p^,^ {b). (2.8) 



From Eqs. (2.3 2.4) one obtains : 



O\ = -\o\ = \[-Gmx{0)- Ge2{0) + Gm3 (0) + ^'a} . 



(2.9) 



We may note that for a spin-3/2 particle without internal structure, for which Gmi(0) = 3eA, 
Ge2{0) = —3eA, and Gm3(0) = —e^ the electric octupole moment of the transverse charge 
densities vanishes. 



3. Lattice Techniques 



The matrix element of the electromagnetic current, , between two A-states can be decom- 
posed in terms of four independent covariant vertex function coefficients, ai{q^), a2{q^), ci{q^) 
and C2{q^), which depend only on the momentum transfer squared = {pf — piY ijl^: 



{^{Pf,sf)\v^\^{Pi,si)) 



Ua{pf,Sf)ff'^^^Ur{pi,Si) 



(3.1) 



air 



Im-A 



+ 



4ml 



2mA 



Ea and tua denote the energy and the mass of the particle, pi (pf) and Sj (sf) are the initial (final) 
four-momentum and spin-projection, while P = pf + pi. Every vector-component of the Rarita- 
Schwinger spinor U(,{p,s) satisfies the free Dirac equation. Furthermore, two auxiliary conditions 
are obeyed: y^ju^ {p,s) = and P(jU^{p,s) = 0. The vertex function coefficients are linked to the 
phenomenologically more interesting multipole form factors Geo, Ge2, Gmi and Gm3 by a Unear 
relation [|l0|]. The dominant form factors are the electric charge, G^o. and the magnetic dipole, 
Gmi, form factors. 

The interpolating field Xaai^) = ■j^e"''' 2 {u"^ {x)Crad''{x)) u^(;c) + {u"^ {x)Crau''{x)) d^(x) 
has the A+ quantum numbers where C is the charge conjugation matrix. To facilitate ground- 
state dominance a covariant Gaussian smearing on the quark-fields entering Xaa i^) used: 
qp{t,x) = ^^[l + aH{x,y;U)]" qp{t,y),H{x,y;U) =1^1^, ([/^(x,0%-A + t^^lC^" A,0%+a) ■ 
Here q is the local quark field (i.e. either u or d), q is the smeared quark field and is the SU (3)- 
gauge field. For the lattice spacing and pion masses considered in this work, the values a = 4.0 
and n = 50 ensure ground state dominance with the shortest time evolution that could be achieved. 

3.1 Correlation functions 

We specialize to a kinematical setup where the final A-state is at rest {pf = 0) and measure the 
two-point and three-point functions 



Ga,{T\p,tf-U) = Y^e-'^'fPTl,^{x^a{tf,Xf)x.aitiM 



(3.2) 
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X,Xf 



(3.3) 



where is the electromagnetic current, which for Wilson fermion is taken to be the symmetrized, 
lattice conserved current. We work with a representation of the Clifford-algebra in which 74 is 



diagonal. In this representation T'^ - 
o^^^ being the Pauli matrices. The ratio 




and Y^ = \ 




, with k = 1 , . . . , 3 and 



Ga\ir,q,t) Gkk{r^Pi,tf-t)Gkk{r\o,t)Gkk{r\o,tf) 
Gkk{r\o,tf) V Gkk{r\o,tf-t)Gkk{r\pi,t)Gkk{r\pi,tf) 



(3.4) 



with implicit summations over the indices k with k = I, . . . ,3, becomes time independent for large 
Euclidean time separations tf — t and t —ti: 



R„\ir,q,t)^u^\ir,q)^ 



Ota 



(3.5) 



The traces act in spinor-space and the Euclidean Schwinger-Rarita spin sum is given by 



Actt(p) 



2mA 



3 3mi 



(3.6) 



Since we are evaluating the correlator of Eq. ( |3.3| ) using sequential inversions through the 
sink [12 1, a separate set of inversions is necessary for every choice of vector and Dirac-indices. 
The total of 256 combinations is beyond our computational resources, and hence we concentrate 
on a few carefully chosen combinations given by 



i.k.l=i k=\ jM,l=l 

(3.7) 

From these all the multipole form factors can be optimally extracted. For instance the first relation 
in Eq. (3^) is proportional to Gmi, while the third isolates Gei for jx =A. 

3.2 Data analysis 

For a given value of q^ the combinations given in Eqs. (3.7) are evaluated for all different 
directions of q resulting in the same q^, as well as for all four directions }JL of the current. This 
leads to an over-constrained linear system of equations, which is then solved in the least-squares 
sense yielding estimates of Geo, Ge2, Gmi and Gm3- The estimates are embedded into a jackknife 
binning procedure, thus providing statistical errors for the form factors that take all correlation and 
autocorrelation effects into account. 

The details of the simulations are summarized in Table [l|. In each case, the separation between 
the final and initial time is — > 1 fm and Gaussian smearing is applied to both source and sink 
to suppress contamination from higher states having the quantum numbers of the A(1232). For the 
mixed-action calculation, the domain-wall valence quark mass was chosen to reproduce the lightest 
pion mass obtained using Ne = 2 + 1 improved staggered quarks [13, p^]. 
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Table 1: Lattice parameters and results. A^conf denotes the number of lattice configurations, (r^) gives the 
charge radius, is the A+ magnetic moment in nuclear magnetons and is the A+ quadrupole moment. 

2 



^conf 


[GeV] mA [GeV] y/{r^) [fm] 




7 




Quenched Wilson, 32^ x64, a = 


0.092 fm 




200 
200 
200 


0.563(4) 1.470(15) 0.6147(66) 
0.490(4) 1.425(16) 0.6329(76) 
0.411(4) 1.382(19) 0.6516(87) 


1.720(42) 
1.763(51) 
1.811(69) 


0.96(12) 
0.91(15) 
0.83(21) 


Nf = 2 Wilson, 24^ x 40(32 for lightest pion), a = 0.077 fm 


185 
157 
200 


0.691(8) 1.687(15) 0.5279(61) 
0.509(8) 1.559(19) 0.594(10) 
0.384(8) 1.395(18) 0.611(17) 


1.462(45) 
1.642(81) 
1.58(11) 


0.80(21) 
0.41(45) 
0.46(35) 




Nf = 2 + 1, Mixed action, 28^ x 64, a 


= 0.124 fm [15] 



300 0.353(2) 1.533(27) 0.641(22) 1.91(16) 0.74(68) 



4. Results 

In Figs ^ § and ^ we show lattice results for Nf = and Nf =2 Wilson fermions and in the 
mixed action. For the pion masses considered in this work there is agreement among results using 
different actions, with statistical errors being smallest in the quenched theory. 
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S 0.6 

a 






i 2 
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0.4 






1 


0.2 












1 1.5 
in GeV^ 




□ = 691 MeV, dynamical Wilson. 
^ m^ = 509 MeV, dynamical Wilson- 
Q m = 384 MeV, dynamical Wilson! 
m^ = 353 MeV, hybrid 




563 MeV, quenched Wilson. 
490 MeV, quenched Wilson 
41 1 MeV, quenched Wilson' 



CHop 



0.5 1 1.5 

in GeV^ 



Figure 1: Comparison of the three different QCD 
lattice calculations for the A+(1232) form factor 
Geo- The upper curve shows a dipole fit to results 
using the mixed action. Results are also shown for 
two dynamical Wilson fermions. The lower curve 
is a dipole fit to results in the quenched approxi- 
mation []l6|]. The error bands are calculated using 
a jackknife analysis of the fit parameters 



Figure 2: Comparison of the three different QCD 
lattice calculations for the A+(1232) magnetic 
dipole form factors Gm\- The lines show the fits 
to an exponential form. The rest of the notation is 
the same as that in Fig. 1 . 
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691 MeV, dynamical Wilson- 
509 MeV, dynamical Wilson 
384 MeV, dynamical Wilson 
353 MeV, hybrid 



□ = 563 MeV, quenched Wilson 
= 490 MeV, quenched Wilson 
Q m^ = 41 1 MeV, quenched Wilson 



0.5 1 
in GeV^ 



Figure 3: Comparison of the three different 
QCD lattice calculations for the A+(1232) elec- 
tric quadrupole form factors Ge2- The notation is 
the same as that in Fig. 1 . 
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* Hybrid 












1 




1 1 





0.0 



0.1 



0.2 

m 



0.3 

' (GeV^) 



0.4 



0.5 



Figure 4: The magnetic dipole moment in nuclear 
magnetons. The value at the physical pion mass 
(filled square) is shown with statistical and sys- 
tematic errors [^. The solid and dashed curves 
show the results of ChEFT with the theoretical er- 
ror estimate [ p^ . 



P^/2^ Pmi' Ulfm^] 




Figure 5: Quark transverse charge densities in a A+(1232) which is polarized along the positive x-axis. Left 
panel : p^3 . Middle panel : i . The light (dark) regions correspond with largest (smallest) values of 
the density. To better see the deformation, a circle of radius 0.5 fm is drawn for comparison. Right panel 
compares the density along the y-axis for p\ (dashed curve) and [ (solid curve). For the A e.m. FFs, the 
quenched lattice QCD results are used. 



The magnetic moment as function of is shown in Fig. 0, together with a comparison to 
a chiral effective field theory (ChEFT) result |[T7|]. The ChEFT result has one free parameter (a 
low-energy constant) that has been fitted to lattice data, shown by the central line. We also estimate 
the uncertainty of the ChEFT expansion (expansion in pion mass) by the error band in Fig. Q The 
uncertainty of the ChEFT calculation vanishes in the chiral limit because in this limit one simply 
has the value of the low-energy constant (LEC), which lies within the broad experimental error 
band jUa+ = 2J^\ '^{stat.) it \.5{syst.) ±3.0{theory)n!^ In this work, we do not consider the 
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uncertainty in the fit value of tlie LEC due to tlie lattice errors, as the calculations are still performed 
for pion masses where the A is stable (on the right side of the kink). A calculation for pion mass 
values where the A becomes unstable will be a challenge for future calculations. The A moments 
using an approach similar to ours are calculated only in the quenched approximation [18, ^ ]. 
Our magnetic moment results agree with recent background field calculations using dynamical 



improved Wilson fermions plQ. The spatial length of our lattices satisfies L^mTi > 4 in all cases 
except at the lightest pion mass with Nf = 2 Wilson fermions, for which Lstrij^ = 3.6. For that 
point, the magnetic moment falls slightly below the error band, consistent with the fact shown in 



Ref. [21] that finite volume effects decrease the magnetic moment. 




.6 





-1.5 -I -0.5 0.5 1 1.5 
b, Ifm] 



-1.5 -1 -0.5 0.5 1 1.5 
b, Ifm] 




-1.5 -1 -0.5 0.5 
A, Ifm] 



Figure 6: Different field patterns in the quark transverse charge density p^3 in a A+(1232) which is polar- 
ized along the positive jc-axis. Left panel : dipole field pattern; Middle panel : quadrupole field pattern; right 
panel : octupole field pattern. For the A e.m. FFs, the quenched lattice QCD results are used. 



In Fig. |5[ the transverse densities of Eqs. (p3,2.4) are compared for a A+ which has a trans- 
verse spin. It is seen that the quark charge density in a A+ in a state of transverse spin projection 
s± = +3 /2 is elongated along the axis of the spin (prolate deformation) whereas in a state of trans- 
verse spin projection = +1/2 it is elongated along the axis perpendicular to the spin. In Fig. ^ 
the dipole, quadrupole and octupole field patterns for the A transverse densities are clearly seen. 



5. Conclusions 

We have presented a study of the electromagnetic properties of the A(1232)-resonance using 
lattice QCD. The lattice results for the A electromagnetic form factors have been shown here down 
to approximately 350 MeV for three cases: quenched QCD, two flavors of dynamical Wilson 
quarks, and three flavors of quarks described by a mixed action combining domain wall valence 
quarks and staggered sea quarks. 

We have also worked out the formalism which allows to interpret these results in terms of 
the quark charge densities. More specifically, we have established the relation between the light- 
front helicity amplitudes and form factors for the case of electromagnetic interaction of a spin-3/2 
particle, which allowed us to obtain the quark transverse charge densities of the A using lattice 
results for the form factors. 
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The light-front formaUsm allows for a consistent relation between the internal structure of the 
particle and its shape. Our lattice results show that the quark charge density in the A"*" is elongated 

along the axis of the spin. 
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